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We propose to use generic Chern numbers for a characterization of topological insulators. It is 
suitable for a numerical characterization of low dimensional quantum liquids where strong quantum 
fluctuations prevent from developing conventional orders. By twisting parameters of boundary 
conditions, the non-Abelian Chern number are defined for a few low lying states near the ground 
state in a finite system, which is a ground state multiplet with a possible (topological) degeneracy. 
We define the system as a topological insulator when energies of the multiplet are well separated 
from the above. Translational invariant twists up to a unitary equivalence are crutial to pick up only 
bulk properties without edge states. As a simple example, the setup is applied for a two-dimensional 
XXZ-spin system with an ising anisotropy where the ground state multiplet is composed of doubly 
almost degenerate states. It gives a vanishing Chern number due to a symmetry. Also Chern 
numbers for the generic fractional quantum Hall states are discussed shortly. 



A crucial role of phases is one of intrinsic features of 
quantum mechanics and has a long history of investiga- 
tion. Among them, those which have intrinsic geomet- 
ric origins are now understood as geometrical phases 0. 
Aharonov-Bohm effects and Dirac monopoles are typical 
and classic examples where the geometrical phases are 
fundamental. Geometrical features of gauge theories are 
another prototype 0- Also a discovery of Berry's phases 
reveals that the geometric phases and the gauge structure 
are closely related and derived by restricting a physical 
Hilbert space|l|i|. 

On the other hand, with an idea of order parameters, 
symmetry breaking is one of the most fundamental con- 
cepts in modern physics. Quite successfully, this stan- 
dard setup can characterize most of ordered states and 
describe phase transitions and critical phenomena. How- 
ever in low dimensional quantum systems, such as elec- 
trons with strong correlation and spins, quantum fluctu- 
ations prevent from developing conventional orders even 
at a zero temperature. In these systems, quantum phases 
of manybody ground states vary wildly in space and time, 
which destroy the standard orders. 

Of course, the wild quantum phases are not random 
but obey some hidden restriction rules and reflect fea- 
tures of the quantum mechanical wave function. Some of 
them are well known today such as Marshall sign rules 
in the spin systems ^ and fractional statistics of quasi- 
particle (hole) wavefunctions in the fractional quantum 
Hall effect (FQHHGl. String order parameters in the Hal- 
danc spin chains[3BlPl ^-nd the quasi off-diagonal order 
m the FQHjj are also discussed based on the feature. 
Generically quantum states with the characteristic geo- 
metric phases are considered to possess non-trivial topo- 
logical orders [llj. 

Recently wide variety of interesting and physically 
important phenomena have been understood based 
on a concept of the topological order. Some of 



them include quantum Hall effects Iisj . solitons in 
polyacetylens(lJ|, adiabatic transports of charge and 
spins 1 151, itinerant magnetism and spintronics|l(ij and 
anomalous Hall conductances |l7lll8lll9l |. polarizations in 
insulators |23|. the two dimensional carbon sheets [2lll2^ . 
anisotropic su per conductors [2^ E^ . and string-net 
condensations |25j. They are under active studies. 

A local phase of the manybody wave function is ar- 
bitrary but there is some correlation with the phases of 
its neighbors, which brings some gauge structures. In 
these view points, templates of such systems are quan- 
tum Hall states, especially integer quantum Hall (IQH) 
states. There are apparently different QH states with dif- 
ferent quantized Hall conductances. However, any sym- 
metries are not broken among the states but they are 
clearly different physical systems. These states are char- 
acterized by the quantized Hall conductances which have 
an intrinsic topological origin [l^. The topological origin 
of the Hall conductance is clear by the Chern number 
expression [1^ EH . Based on the observation, we pro- 
pose to characterize the topological orders by the generic 
Chern numbers\2^. The Chern numbers in topological 
ordered systems are kinds of order parameters in conven- 
tional ordered phases. In the same manner as the usual 
phase transition is characterized by a sudden change of 
order parameters, topological phase transitions are char- 
acterized by a discontinuous change of the Chern num- 
bers. 

To have a well defined Chern number, we need an ex- 
istence of a generic gapH. Topological insulators are 
defined as physical systems with this generic energy gap. 
Then the Chern numbers are always integers and the 
topological phase transition is characterized by a discon- 
tinuous change of the Chern numbers which are always 
integers. This integral property of the Chern number 
implies a stability of the characterization against a small 
perturbation. However, as in the case of the edge states of 
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the quantum Hall effect [29| and Kennedy's triplet states 
in the Haldane spin chains |30|. the topological ordered 
state are quite sensitive to a geometrical change of the 
physical system such as an existence of edges and bound- 
aries. It contrasts with the conventional order where 
boundary conditions are always negligible in the ther- 
modynamic limit. Therefore a translational invariance 
is fundamentally important to describe topological or- 
dered states. On the other hand, in many cases, as far as 
physical observables are concerned, a topological order is 
hidden in a bulk and only reveals its physical significance 
near boundaries of the system. 

Generically speaking, to define the Chcrn numbers C 
for a physical (many particle) wave function, ?/;, we need 
to require the wave function to depend on multiple pa- 
rameters, X E V, dimV > 2. Most common such pa- 
rameters without disturbing bulk properties are multiple 
Aharonov-Bohm fluxes on a genus g Riemann surface. 
When the topological order is non-trivial, there can be 
inevitable topological degeneracies 111, such as a g^-fold 
degene racy of the FQH state with a filling factor 1/g on 
a torus 0[ 1^. The degeneracy of a generic ground 
state will be discussed later. Here we just point out that 
one has to consider non-Abelian gauge structures arising 
from it 0, . This is crucial for a numerical concrete 
characterization of the topological insulators. Especially 
an explicit gauge fixing for the degenerate multiplet is 
required to perform calculations. j28j. 

Quantum Spin Systems as Topologi- 
cal Insulators : To describe a characteriza- 
tion of the topological order, let us consider a 
generic translational invariant spin- 1/2 hamilto- 
nian on a d-dimcnsion orthogonal lattice — 

r^m ^ ymi...ym,^ ^ ^ (mi,---,mrf) and 
*S{r) = XS'-'{r),Sy{r),S''ir)) is a spin-1/2 opera- 
tor at a lattice site r and h( is a local hamiltonian 
which depends on several spins at ri, r2, ■ ■ ■ . It gener- 
ically breaks several symmetries explicitly such as a 
parity, a chiral symmetry, and a time reversal symme- 
try. The operator is a translation in /i-direction, 
r^S'(r)T^ — S{r + a^^) (a^ is a unit translation in 
jji direction). We use a periodic boundary condition 
T^f^ = 1 (to^ — 0, • • ■ , d) to avoid disturbing bulk 
properties by possible edge states. We propose to use 
twisted boundary conditions for the spin model and take 
the twists as the parameters x as discussed below. 

Local Gauge Transformation and Twists : Let 
us consider a local gauge transformation of a string type, 
that is, local spin rotations at a unit cell label m as 
S^(r^) = Q(7)5(r^) with 3x3 matrix Q(7) = e^^, 
7 = m-6> where X"'^ = im^Tr cr"cr''cr^, 6» (6*1, • • • , 61^) 
and n = {nx,ny,nz) (|n| = 1) is a fixed rotation axis. 
Also 77 is a label to distinguish intra unit cell spins. 
The simplest example is given by taking n = (0,0,1) 



(cos 7 — sin 7 \ 
sin 7 cos 7 with "f — m ■ 0. We 

further assume the local hamiltonian hi is of the gauge 
interaction type as 

/i,(S(ri), 5(r2), • . • ) =/i«(0 = 0; {r, - r,}; S{r,), S^r^), ■■■) 
=h''{e-{r,^r,}-S's{ri),S',{r2),---) 
^h%S'o{r,),S'oir2),---) 

with some function h'-^. That is, the twisting parameters 
only affect the hamiltonian through the relative positions 
of the local spins. Examples of such interactions for the 
above rotation around z-axis are h^^" = *S^(ri) JS(r2) 
with J = Jdiag(l,l,A) and hf = JcS{ri) ■ {S{r2) x 
Sirs)) = Jce^JkS,irl)SJ{r2)Skir3)^. They trans- 
form respectively as h^'"''' = ^{e-'^-'^^-'^^'' Sij^S2g^ + 
h.c.) + XS4^S2l, /i^^ = %5i^^(ze-(«-«3)52+^3f^ + 
h.c.) + (cyclic perm.), 6i — rrii ■ 9, i — 1,2,3 where rrii is 
a unit cell labeling of the spin S{ri). 

The hamiltonian is periodic in the original rep- 
resentation by S^s but is not periodic in the one by 
twisted Se's as S'iT^'-r) = T'l''S'{r) with T'^" = 
exp(-n^6'^L^X). 

Now let us define a translational invariant twisted 
hamiltonians ff'^ by a representation by the twisted S'g 

as H^{e) = Y.mT''^Ki.S'e{r,),S's{r2),---)r'^^w\th 
a periodic boundary condition T'^^ = 1. In the original 
spin operators S, is given by with the twisted 
boundary condition = exp(-|-0^L^JC). 

The two hamiltonians, and , are both transla- 
tionally invariant in representations by S and S'g respec- 
tively. One may expect an macroscopic 0{V) energy 
difference between their ground state energies. How- 
ever, as discussed, the contribution should be at most 
C(|9p^|) due to the gauge invariance where is a 

(hyper) area of the system perpendicular to the r^^-axis 
where V — Li ■ ■ ■ Ld- That is, the difference of the en- 
ergy should be a finite size effect. Thus the difference 
between and is negligible in the thermodynamic 
limit y — > 00 when we discuss the bulk properties in a 
usual manner. 

Another important point for the present construction 
is that the twisted hamiltonian is translational in- 
variant in the So representation. Then edge states never 
appear in any representation, which is especially impor- 
tant to pick up only bulk properties through probes by 
twisted boundary conditions. 

Degeneracies and Ground State Multiplet : A 
topological order on a non-zero genus Riemann surface is 
one of the reasons for the ground state degeneracy, which 
is the topological degeneracy The simplest example 
is just a manybody state with two-dimensional periodic 
boundary conditions!^ 113- Also if the system has a 
standard symmetry breaking, such as ising orders, a fi- 
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FIG. 1: Schematic spectral flow with parameters x. 



nite system has almost degenerate ground states corre- 
sponding to linear combinations of the symmetry broken 
states 35]. For a finite system, the degeneracy can be 
lifted and the splitting is estimated as « e"*^^ for the 
symmetry broken states. Also if the ground state has a 
finite spin moment which may not be macro as a ferro- 
magnet, there occurs a spin degeneracy. Some of these 
degeneracies can be approximate for a finite system and 
may be sensitive to the boundary condition and twisting 
parameters, such as 0^'s. In these cases, the lowest en- 
ergy gap is not a physical one and may vanish in the ther- 
modynamic limit. The physical energy gap for the bulk is 
an energy gap above these almost degenerate states. We 
define a ground state multiplet by a collection of these al- 
most degenerate states near the ground state and define 
a Chern number for this ground state multiplet. (See the 
FigP) (361] Since the two hamiltonians and H^{9) 
differ only boundary terms, bulk properties of the two 
should be the same. Then, for the topological insulators, 
the energy gap above the ground state multiplet is stable 
against perturbations. If the ground state multiplet is 
well separated from the above in a finite system, we do 
not need to take the thermodynamic limit. 

The twist we proposing is a boundary perturbation 
in a particular representation. However it also preserve 
the translational symmetry up to a unitary equivalence. 
Then, based on a discussion of edge state picture, we ex- 
pect an energy separation of the topological degeneracy, 
that is, a band width of the multiplet by the twist as 
e~^/4 where L is a minimum linear dimension of the sys- 
tem and ^ is a typical length scale of the ground state 
multiplet. It can be different from the conventional bro- 
ken symmetric cases. 

Chern Numbers for the Spins : Let us de- 
fine a total parameter space by V = {(6*1, • • ■ ,6d)\0f_i £ 
[0, 27r/ip]}. Since exp(27rX) = J3, we have 
H^{e + (•••, 0, + 27r/L^, 0, • • • ) = in the S- 

representation. Then the twisted hamiltonian H^[0) 
is well defined on V without boundaries as V = T*^. 
I37I ] Any two dimensional integration surface S{c V) 
without boundaries is used to define the Chern num- 
bers. S = Tfj = is the simplest example. 



Now define a ground state multiplet ^{x),x £ V. It 
is a iV X q matrix as '^{x) = (■(/'i(a;), ■ • ■ , '0ij(2^))with 
H^{x)'il)j{x) = ej^j{x),j = 1,2,---, < ej,(i < j), 
where V'j is a column vector in a many spin Hilbert space 
with a dimension N and g is a dimension of the ground 
state multiplet. The generic energy gap condition for the 
multiplet is given as eq(x) < eq+i(a;),Vx G S. This is a 
definition of the topological insulators. 

Define a non-Abelian connection one-form A which 
is an g X g matrix as ^ = d'^ and a field strength 
two- form ^ ~ dA. + A?. The first Chern number [2] is 
then defined by = 2^/sTr:F = ^J^TidA. It 
is a topological integer which is stable against pertur- 
bation unless the generic gap collapses. We use these 
integers depending on a choice of S to characterize the 
topological orders. Changing a basis within the mul- 
tiplet space, '^''(a;) = '^{x)u;{x), {ljuj^ — Iq) gives a 
gauge transformation A! = ^'d^' = uj^^Auj + u)^^du3 
and :F' = u3-^J^u3 H IH]. The Chern number is a 
gauge invariant but we need to fix the gauge to evalu- 
ate this expression [2^. Take a generic arbitrary multi- 
plet, and define an overlap matrix as 0$ = $^P$ 
where P = '4"4't is a projection into the ground state 
multiplet which is a gauge invariant. Then define re- 
gions 5*, R — 1, 2, ■ • • , as (infinitesimally) small neigh- 
borhoods of zeros of det 0$ {x) and <S* as a rest of 
S. Then the first Chern number is written as Cs — 

The field n is defined as = Arg det*^P* = 

Arg det rj — Arg det r) where ^ is also another generic ar- 
bitrary multiplet, r] = "Jf^^ and f) = '5'^$. The matrices 
Tj and fj depend on the choice of the multiplet ^ but the 
difference of the arguments is a gauge invariant. 

The field depends on a choice of $ and ^ but the 
total vorticity Nq{S) is a gauge invariant and indepen- 
dent of the choice. The field fl reflects a phase sensitivity 
of the multiplet by the twist when one fixes $ and It 
is illustrative to show and it supplies information of 
the ground state multiplet. Also when the integration 
surface S is contractible to a point keeping a generic en- 
ergy gap, the Chern number vanishes from a topological 
stability. 

Ex.l:Two-Dimensional Spin Model : The present 
formulation can be effective for characterization of topo- 
logical ordered phases in any dimensions, such as chi- 
ral spin states l34|. To have a finite Chern number, one 
needs to break time reversal symmetry as for the quan- 
tum Hall states [s^. The simplest example of hi can 
be a sum of local pair-spin interactions /i^**" with a 
symmetry breaking term hf' discussed above as hi = 

J2pair^e^" + ^f'- Here let us just show an example 
with a degeneracy to show the present general procedure. 
Considering only a nearest neighbor exchange interac- 
tion and assume n — (0,0,1), the local hamiltonian is 
given as /i« = J{'^ie-^'' SglSg^+e^'" Se''^S9+) + 



4 




(b) 2II/L, 



vorticity of at the zeros of | detT?|^pq. 
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FIG. 2: (a)Three lowest energies of the XXZ model on a 
4x4 square lattice with twists in the total Sz = sector is 
shown. (A = 1.3) (b) A field Q, of the ground state multiplet 
composed of the lowest two eigen states for some choice of 



XSe'iSg,) where = SgiT^r) and Sg ^ Se{r). In 
the case, the twisted boundary condition for H^{0) in 
S is given by the following matrix — Q^(7) with 
7 = OxLx + OyLy. Then wc can take a 2 dimensional 
torus T2 = {{e^,ey)\Q < Of, < 27r/L^} for the integra- 
tion surface S. This is a nearest neighbor XXZ model 
on a square lattice with twists. With an ising anisotropy, 
A > 1, the ground state of an infinite system has a long 
range order and has a finite energy gap. We show nu- 
merical results for a system with J — 1 and A = 1.3. 
The ground state of a finite size system is given by a 
bonding state between two symmetry broken states with 
antiferromagnetic (ising) order. The next lowest state 
is an anti-bonding state of them and the energy separa- 
tion between is expected to be oc e"^^^*/^ where f is 
a typical length scale. A physical ising gap to flip one 
ordered spin is given the one above, that is, the second 
lowest one. Therefore the ground state multiplet is com- 
posed of the two low lying states including the finite size 
ground state. (See Fig[2Ia)) The field fl of the two di- 
mensional ground state multiplet is shown in Figl^Jb). 
As discussed in the reference |28f. the Chern number is a 
sum of the vorticity at the zeros of detO$ = | detrjp. 
In the present example, it is as expected for a chiral 
symmetric system. 

Ex.2, Manyparticle States in the First 
Quantized Form : The same procedure is 
also applied for a manyparticle state in the first 
quantized form, such as the generic FQH States 
\E'fe(x; ri, cTi; • • ■ ; r^, <tjv) where k denotes a label of the 
(topological) degeneracy of the ground state multiplet 
and X = {9x,0y) is a set of parameters specifying 
twisted boundary conditions on a torus 32] . By taking 
a reference multiplet as $^(ri, cri; • • • ; tat, cttv) — 
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the degenerate multi- 
the field Q,{x) = Argdet 7)^77 
*fe(a;; , cr?; • • ■ ; rl, crl) and 
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IS given 

{^Ijfc = '^k{x]r{,(T{\--- ;r^,(T^) where fc, ^, ^ run over 
{1, • • ■ ,9}. The Chern number is evaluated as a total 
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